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We present new mass models for the gravitational lens system B 1938+666, using multi- 



& \ wavelength data acquired from Keck adaptive optics (AO) and Hubble Space Telescope (HST) 

observations. These models are the first results from the Strong-lensing at High Angular Res- 
olution Program (SHARP), a project designed to study known quadruple-image and Einstein 
ring lenses using high-resolution imaging, in order to probe their mass distributions in un- 
precedented detail. Here, we specifically highlight differences between AO- and ffiT-derived 
lens models, finding that - at least when the lens and source galaxies are both bright and 
red, and the system has a high degree of circular symmetry - AO -derived models place sig- 
nificantly tighter constraints on model parameters. Using this improved precision, we infer 
important physical properties about the B 1938+666 system, including the mass density slope 
. of the lensing galaxy (7 = 2.045), the projected dark matter mass fraction within the Einstein 

radius (Mdark/^iens = 0.55), and the total magnification factor of the source galaxy (~ 13). 
Additionally, we measure an upper-limit constraint on luminous substructure (My > 16.2), 
based on the non-detection of bright satellite galaxies in all data sets. Finally, we utilize the 
improved image resolution of the AO data to reveal the presence of faint arcs outside of the 
primary Einstein ring. The positions and orientations of these arcs raise the intriguing possi- 
bility that B 1938+666 has a second source galaxy, located at a more distant redshift. However, 
future work is needed to verify this hypothesis. 
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1 INTRODUCTION 



An understanding of the nature and distribution of matter 
on small (< 1 Mpc) scales is essential to modern astro- 
physics. Measuring the shapes o f galaxy mass profiles reveals 
the presence of dark matter (e.g., iRubin. Peterson, & Ford Il980l: 
van Albada & Sancisil [l986l: Ide Blok & McGauehl 1 19971: iBosmal 



1999: lGavazzi et all 120071: iConroy et al.ll2007l : iDutton et al.ll2011l : 
Ruff et alj|201ll : ISuvu et al.ll2012h and its in t eractions with bary- 
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constrain cosmological parameters (e.g. the baryon fraction Qb) 
and plac e estimates on the efficiency of star-formation in galax- 
ies (e.g., Fukugita. Hogan. & Peebles! [l9 9 8: Hev mans et al J 1 20061 : 
iNapolitano. Romanowskv & TortoraluOlOl : lLagattuta et al.l2010h . 
Observing changes in any of these quantities over cosmo- 
logical ti me is a key component in studying galaxy evolu- 
tion (e.g.,|M andelbaum et al. 2006; Behroozi, Conrov & Wechsleq 



120101 : lLagattuta et al.ll2010h 



While there are many techniques capable of measuring mass 
on small scales, gravitational lensing stands out as an especially 
powerful choice. Unlike other methods, lensing directly measures 
a total (baryonic + dark matter) mass without requiring this mass 
to be luminous or in any specific dynamical state. Furthermore, a 
lensing analysis is not limited to the local Universe, but rather can 
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be applied to systems located over a wide range of cosmological 
distances. Strong gravitational lensing, in particular, ca n provide a 
wealt h of information about the nature of galaxies (e.g., Kochanek 
2006 and references therein). With typical image separations of ~ 
1 arcsec, galaxy-scale strong lenses produce mass estimates close 
to the centre of galaxies, giving information about both the baryon- 
dominated luminous core and the inner regions of the dark matter 
halo. This relatively small angular size (corresponding to physi- 
cal scales of 5-10 kpc at typical lensing-galaxy redshifts) decreases 
the probability of foreground interlopers contaminating the line of 
sight - a problem that can bias mass estimates obtained from group 
and cluster-scale strong lenses, and can strongly dilute the signal 
measured from weak lensing - suggesting that the observed lens- 
ing signal will be dominated by the mass of the lensing galaxy. 

In this paper, we investigate the galaxy- scale gravitational 
lens B 1938+666. First discovered as part o f the Jodrell Bank- 
Very Large Array A s trometric Survey (JVA S; IPatnaik etal 1 ll992l : 
iBrowne et all Il998l : IWilkinson etal] Il998l) . initial radio obser- 
vations of B 1938+666 showed a quadruply-imaged background 
source, configured into a par tial Einstein ring, along with a second , 
doubly-imaged component dPatnaik et a l. 1992; King et al]fl997h . 
Follow-up imaging in the near-infrared (NIR) and optical uncov- 
ered a bright red object that was thought to be the galaxy lens- 
ing the radio emission dRhoads. Malhotra & Kun dic 1996). This 
emission was later shown by NIR Hubble Space Telescope (HST) 
imaging to be composed of the light from both the lensing galaxy 
and a complete E instein ring of the backg round source galaxy 
dKing et alJll998h . lTonrv & Kochanekl l2000l) measured a redshift 
of z\ — 0.881 for th e lensing galaxy from optical spectroscopy, 
while Riechers (2011) determined a redshift of z s = 2.059 for the 
source from CO observations. 

As B 1938+666 has an Einstei n ring, its lensing mass model 
can be determined to high precision. Ko chanek. Keeton & McLeodl 
(l200ll) showed that, when compared with two-image quasar or arc 
lens systems, constraints from an Einstein ring can be used to break 
degeneracies between the monopole moment of the gravitational 
potential and higher order terms. Removing these degeneracies al- 
lows for a more robust measurement of the slope of the lensing 
galaxy's mass profile and, when combined with time delay infor- 
mation, can p rovide an unambi guous measurement of the Hubble 
constant (e.g. ISuvu et aDl2010h . At the same time, with so many 
constraints placed on the model, Einstein ring lenses are sensi- 
tive to pertu rbations from the smooth gravitational potentia l of the 
lens galaxy (lKoopmansll2005l : IVegetti & Koopman s)l2009ah . Satel- 
lite galaxies and dark matter sub-haloes that orbit the main lensing 
galaxy can give rise to these perturbations, and while these objects 
are often too faint to be directly observed, gravitational lensing can 
be used to detect them indirectly through their mass signatures. In 
this way, careful scrutiny of Einstein ring lenses (or those with ex- 
tended gravitational arcs) provides one of the best opportunities to 
reveal the presen ce of, and constra i n the properties of, extragalac- 
tic substructure (|More et all 120091; fVegetti Czoske. & Koopmansl 
l2010l : IVegetti et al.l2010ll2012l : ISuvu & Halkolall2010h . This would 
serve as a direct test of the ACDM numerical simulations that 
predi ct this substructure (e.g.. iDiemand et al"1l2008l; ISpringel et al.l 
l2008h . and would thus advance our understanding into the nature 
of galaxy formation and evolution. 

For the reasons described above, we choose to study 
B 1938+666 using high-resolution, NIR adaptive optics (AO) as 
part of the Strong lensing at High Angular Resolution Program 
(SHARP). This new project aims to obtain high resolution images 
of lens systems to study the mass distributions of lensing galaxies 



between 0.3 < z < 1 to an unprecedented detail with state-of- 
the-art imaging and lens modelling techniques. In this first paper, 
we focus on a comparison between imaging data taken with space- 
and ground-based observatories, in particular, we determine the rel- 
ative benefits for lens modelling. We also present a new smooth 
mass model for our test system B 1938+666 that is used to deter- 
mine the properties of the lens galaxy and give a robust estimate of 
the magnification of the NIR component of the background galaxy. 

This paper is organized as follows. In Section |2] we describe 
the wider goals of the SHARP survey since this is the first pa- 
per in the series. In Section [3] we briefly describe the new and 
archival multi- wavelength imaging data for B 1938+666, new IR 
spectroscopy of the lensed source and the techniques used to re- 
duce them. In Section [4J we present the results of lens modelling. 
We discuss our results in Section|5] Finally, we summarize and con- 
clude in Section[6] Measurements describing the relative fraction of 
substructure within the B 1938+666 len sing galaxy are presented in 
a companion paper dVegetti et al.ll2012h . 

Throughout this paper, we assume a cosmological model of 
Ho = 100 h krns" 1 Mpc -1 , Q m = 0.3, and Q A = 0.7. All 
magnitudes presented in this work are AB magnitudes. 



2 SHARP RATIONALE 

The image separations seen in galaxy- scale lenses are typically on 
the order of an arcsecond. Thus, observing these systems with high 
angular resolution instruments can provide unambiguous estimates 
of individual image positions and magnitudes, information that is 
critical when measuring mass distributions with gravitational lens- 
ing. Furthermore, any improvements in angular resolution com- 
pared to traditional ground-based imaging will provide enhanced 
detectability of faint objects with small angular extent, and im- 
proved sensitivity to small-scale perturbations of the surface bright- 
ness of lensed extended emission. At optical and NIR wavelengths, 
there are currently two techniques for obtaining the necessary data: 
space-based imaging with HST and ground-based AO imaging. The 
SHARP survey utilizes both approaches in order to build up a sta- 
tistically significant sample of lens systems for which deep high- 
resolution imaging has been obtained. 

A particular focus of the SHARP survey is to detect and mea- 
sure the mass of substructures associated with the lensing galax- 
ies, without regard to whether the substructures are luminous or 
dark. Numerical simulations of galaxy formation predict a large 
amount of substructure for a galaxy-mass halo, with a mass frac- 
tion of / su b = 5-10 per cent of the total halo mass contained in 
substructures with masses between 4 x 10 6 Mq and 4 x 1O 9 M 
dDiemand et all 2008; Sp ringel et ai1l2008h . Furthermore, the sim- 
ulations converge on a substructure mass function of dN/dm oc 
m a , where a — —1.9 ± 0.1 (e.g. , IDiemand. Kuhlen. & Madaii 
120071) . IVegetti & Koopmansl d2009bh have shown that a Bayesian 
analysis of a sample of lenses that have been surveyed for substruc- 
ture down to some mass threshold can provide meaningful con- 
straints on / su b and the slope of the mass function, a. Given a sam- 
ple of ~30 lens systems, a mass detection threshold of 3 x 10 8 M© 
or better, and a reasonable prior on a, good constraints on / su b are 
obtained. Note that even the non-detection of substructure down 
to the mass limit provides useful information for constraining / su b 
and a, especially as the mass detection limit becomes smaller. 

We are therefore pursuing two complementary methods for 
detecting substructure in lens systems, both of which utilize deep, 
high-resolution imaging. The first is to detect luminous satellites 
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Figure 1. Multi-wavelength imaging of B 1938+666. Top row: HST images taken using the F555W (V-Band), F814W (I-Band), and F160W (H-Band) filters, 
respectively from left to right. Bottom row: Keck-II Telescope LGS AO images showing the H-band (left panel) and K'-band (middle panel) light. A colour 
composite image (right panel) was made by combining the HST V- and I-band images with the LGS AO H- and K'-band images. All images are oriented 
with North pointing up and East pointing to the left. The lensing galaxy and source galaxy are clearly seen in all of the near-infrared bands, whereas only the 
lensing galaxy can be seen in the optical bands. Additionally, a second set of faint arcs can be seen on the eastern side of the K'-band image (see Section 1531 
for further details). 



directly, thus fixing the location of the substructure. The substruc- 
ture can therefore be included in the mass model o f the lens and 
its mass can be determined (e.g., iMore et al ]|200§. This method 
is most effective for lens systems with either four lensed im- 
ages or extended lensed emission. The other approach is to de- 
tect substructures thr ough their gravitati onal effects on extended 
lensed emission (e.g., Ve getti et al.|[2010h . Thus, the SHARP sam- 
ple consists of lens systems with quadruply imaged quasar lenses 
and lenses with partial or complete Eins tein rings. Results on two 
four-image lenses have been reported by iMcKean et al.l (120070 and 
lLagattuta Auger & Fassnachtl j2010h . Here we report the results 
for B 193 8+666, which has an almost complete Einstein ring. A de- 
tailed description of the SHARP sample will be presented by Fass- 
nacht et al. (in preparation). 



3 OBSERVATIONS AND DATA REDUCTION 
3.1 Adaptive optics imaging 

We observed the B 1938+666 system on UT 2010 June 29 and 30 
with the NIRC2 camera on the Keck-II telescope, using the Laser 
Guide Star (LGS) AO system. We used the narrow camera, with a 



Table 1. A summary of the B 1938+666 optical and infrared imaging obser- 
vations. 



Date 


Telescope 


Instrument 


Filter 


^exp (s) 


1997 Aug 13 


HST 


NICMOS/NIC1 


F160W 


10800 


1999 Apr 24 


HST 


WFPC2 


F555W 


2800 


1999 Apr 24 


HST 


WFPC2 


F814W 


3000 


2010 Jun 29 


Keck II 


NIRC2 LGSAO 


K' 


14760 


2010 Jun 30 


Keck II 


NIRC2 LGSAO 


H 


6600 



field-of-view of 10 arcsec x 10 arcsec, and a pixel scale of 0.01 
arc sec. The tip- tilt correction was obtained through simultaneous 
observation of a magnitude R = 15 star at a distance of 18 arcsec 
from the lens. We observed the system using both the H and K ; 
filters. Details of the observations are given in Table [T] 

The data wer e reduced following the method presented in 
lAuger et all (1201 lh . The final reduced images from the AO obser- 
vations are shown in the lower row of Fig. [T] The images have three 
distinct features: a bright, compact lensing galaxy in the centre of 
the image, a high signal-to-noise ratio (SNR) Einstein ring, and a 
set of faint arc structures outside of the ring on the east side. Since 



© 2012 RAS, MNRAS 000, ??-?? 



4 D. J. Lagattuta et al. 



CD 

U 1.0 

I 0.8 
£ 0.6 
^ 0.4 
o 0.2 
"t5 0.0 
£-0.2 

2.0 

x 1.5 
u_ 

§ 1.0 
S 0.5 
0.0 
-0.5 





N5 



atmospheric transmission 




N7 



16000 



18000 20000 22000 
Wavelength (Angstroms) 



24000 



Figure 2. The NIRSPEC spectra of the combined emission from the lens- 
ing galaxy and background source (solid) and the la error spectrum (dot- 
ted). The spectra have been smoothed with a nine-pixel boxcar, with the 
points weighted by the inverse of their variances. The expected positions of 
the [Oin] AA4960, 5007 and [Sn] AA6716 6731 lines at the background 
source redshift of z = 2.059 (Riech ersboilh are marked, while the vertical 
dashed line represents the expected position of Ha. Additionally, the atmo- 
spheric transmission over the entire wavelength range of interest is shown 
in the upper panel. 



we did not observe any photometric standards during the data ac- 
quisition, we calibrate our photometry by scaling the observed flux 
to match with previously reported data. The total (lens + ring) ri- 
band fl ux is matched to the NICMOS data presented bv lKing et alj 
dl998h (F160W = 19.4 ± 0.3) while the K'-band data are matched 
to the NIRC value reported bv lRhoads. Malhotra & Kundid dl996h 
(K' = 19.0 ±0.1). 



3.2 Archival Hubble Space Telescope imaging 

The B 1938+666 system has previously been observed with the 
HST, at MR (GO 7255; PI Jackson) and optical (GO 7495; PI 
Falco) wavelengths. The optical data were obtained as part of the 
Cf A- Arizona Space Telescope Lens Survey (CASTLES) progranfl 
These observations are also summarized in Table Q] The HST data 
were obtai ned from the archiv e and reduced using the routines de- 
scribed in lAuger et all d2009l) . The final images are presented in 
the upper row of Fig. Q] and, in the case of the NIR imaging, show 
the same structure that was observed in the Keck AO data. In the 
F555W and F814W bands, the Einstein ring was not detected. 



3.3 Near infrared spectroscopy 

The B 1938+666 system was observed with the near-infrared 
echelle spectrograph (NIRSPEC: [McLean et al.ll998l) on the Keck- 
II telescope on the night of UT 2006 August 14. The observing con- 
ditions were marginal; scattered clouds throughout the night signif- 
icantly altered the sky transparency, which led to poor data qual- 
ity. The spectra were obtained through the NIRSPEC-5 (roughly 



http://www.cfa.harvard.edu/castles/ 



H-band) and NIRSPEC-7 (roughly K-band) filters. Four exposures 
were taken in each band, with each exposure consisting of three co- 
adds of 300 s. Due to the poor observing conditions, only two of 
the NIRSPEC-5 exposures were usable. All four NIRSPEC-7 ex- 
posures were used. The spectra were reduced with a Python-based 
pipeline, which consists of steps to subtract the sky emission, re- 
ject cosmic rays, rectify the two-dimensional spectra, wavelength 
calibrate, and extract the object spectra. The final spectra have 
pixel scales and spectral resolutions of 3.56 A pix -1 and ^1600 
in NIRSPEC-5 and 4.08 A pix -1 and -1400 in NIRSPEC-7. 

The source redshift was unknown at the time of the observa- 
tions, and both the Ha and H/3 lines fall outside the wavelength 
range covered b y the spectra. H owever, given the source redshift 
of z s = 2.059 kiechersll20TTh . the [Om] AA4960, 5007 emis- 
sion features should fall in the NIRSPEC-5 spectrum and the [S II] 
AA6716, 6731 features should fall in the NIRSPEC-7 spectrum. 
The spectra are shown in Fig. [2] They have been normalized and 
then smoothed by a 9-pixel moving average, with each point be- 
ing inverse-variance weighted. Even with the smoothing, no clear 
features are seen in the spectrum. However, the expected locations 
of the [O III] and [S II] emission lines (marked in Fig. |2} do corre- 
spond to weak peaks in the spectrum. Thus, although we would not 
have been able to unambiguously measure the source redshift with 
the NIRSPEC spectra, w e can say that ou r data are consistent with 
the redshift measured bv lRiechersI (1201 ll) . 



4 THE LENS MODEL 

Strong gravitational lenses with extended- source structures are fre- 
quently modelled by first determining and subtractin g the surface 
brigh tness distribution of the foreground galaxy fe.g. jBolton et alj 
2006), although in cases where the background object is very bright 
this can lead to an over- subtraction, where parts of the strongly 
lensed features are fitted and remove d along with the foreground 
galaxy light (e.g. lAuger et all l201ll) . At NIR wavelengths, the 
(lensed) B 1938+666 source galaxy is roughly equal in brightness 
to the foreground (lensing) galaxy, so the system presents a similar 
modelling challenge. Therefore, we first fit a simply parametrized 
lens model to the data in order to quantify and remove the fore- 
ground galaxy light. We then fit a more detailed lens model to the 
residual data, in order to precisely infer the properties of the mass 
distribution. 



4.1 Surface brightness modelling 

Modelling the surface brightness distribution for this system re- 
quires that we also determine an approximate model for the lensing 
mass distribution, in order to disentangle the foreground lens and 
background source light. To do this, we take an approach similar to 
lAuger et all d201 ll) and employ an elliptical power law mass model 
(e.g. lBarkanall 19981) to describe the mass distribution, while includ- 
ing an external shear contribution. The foreground and background 
galaxy su rface b r ightne ss distributions are modelled with (possibly 
multiple) ISersid dl963h profiles. The models are fitted to the data 
using an adaptive simulated annealing scheme, and each potential 
model includes either a single surface brightness component for 
both the lens and source galaxies, a single surface brightness com- 
ponent for one galaxy and two components for the other, or two 
surface brightness components for each galaxy. 

The foreground galaxy is well-modelled with a single Sersic 
component, while the background source modelling strongly 
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favours two components. Although we find that there is signifi- 
cant covariance between the structural properties (i.e., the Sersic 
indices, effective radii, and total magnitudes) of the foreground and 
background components when fitting to the AO images, the central 
surface brightness distributions are robustly segregated. For exam- 
ple, the total magnitude of the foreground galaxy can change by 0.5 
magnitudes if one or two components are used for the background 
source, but the flux within the Einstein radius - the radius of the 
ring produced when a part of the lensed object sits directly behind 
the lensing galaxy - only varies by about 10 per cent (~0.1 magni- 
tudes); the covariance between the total galaxy properties is likely 
due to inadequate knowledge of the AO PSF. We therefore restrict 
our discussion of the surface brightness properties to the inferred 
luminosity within the Einstein radius and the position of the source 
relative to the lens. The aperture magnitudes within the Einstein 
radius are 21.0 ± 0.15 and 20.6 ± 0.10 in the H- and K' -bands, re- 
spectively. Note that with lensing we measure the total mass prop- 
erties within the Einstein radius and we are therefore able to draw 
robust conclusions about the relationship between mass and light 
in the centre of the lensing galaxy, which we discuss in Section l5^2l 

4.2 Mass Modelling 

We model the Keck AO H- and K'-band data sets along with the 
HST NICMOS F160W-band data set using the adaptive a nd grid- 
based Bayesian technique of Vege tti & Koopmansl l2009al ). An ex- 
tensive description of the modelling procedure in the context of 
Bayesian evidence optimization is provided there. In short, we 
model the system by assuming an elliptical power-law mass dis- 
tribution [p(r) oc r -7 ] for the lens galaxy, which we use to recon- 
struct the surface brightness of the source galaxy (after first sub- 
tracting off the lens galaxy light profile; see Section |4J1 on a two- 
dimensional grid. This grid is adaptive in magnification and built 
by casting pixels from the image plane back to the unlensed source 
plane. Due to the high resolution of the AO data, we only cast one 
point out of each six-by-six pixel sub-grid back to the source plane, 
whereas for the lower resolution HST data, we use a smaller three- 
by-three pixel sub-grid. Note that the full pixel grid is used when 
the reconstructed source surface brightness is cast forward to the 
image plane and the residuals are calculated. 

The best lens models (as determined from maximum- 
likelihood analysis) for the K'-band, H-band, and NICMOS data 
sets are presented in Figs. |3(a)||3(b)| and |3(c)| respectively, while the 
individual parameters of these models are listed in Table In addi- 
tion to the maximum-likelihood values, we also derive mean values 
and confidence intervals for each parameter's marginalized poste- 
rior probability distribution function by exploring the Bayesian evi- 
dence in the full multidimensional parameter space (see Table |2]for 
the results). One- and two-dimensional slices of these marginalized 
posterior p robability distributions a re shown in Fig. @] Using Multi- 
nest v2.7 dFeroz & Hobsonl 120081) . we are able to integrate over 
this posterior probability and obtain the marginalized Bayesian ev- 
idence, which is the probability of the data given the model family. 

Incorporating a robust PSF model can be challenging for AO 
data, since the PSF varies rapidly over time. This is a particular con- 
cern for the B 1938+666 AO data sets as we are unable to observe 
the lens system and a PSF star simultaneously. Therefore, we carry 
out the modelling using different PSF stars taken at different times, 
leading to a different "best" model for each PSF. As the arc is suf- 
ficiently extended and the dynamic range is relatively low, we find 
that the lens modelling o f B 1938+666 is not significantly affected 
by the choice of PSF (see Ve getti et"aDl2012l for details). However, 



the Bayesian evidence allows one to choose the best PSF model 
objectively. For example, the PSF used for the model labelled Mk 
in Table |2]yields the largest evidence value for the K'-band data. 

The best models for each data set are consistent, as can be seen 
by comparing the maximum-likelihood parameters given in Table 
[2 Many of the marginalized posterior probability distributions are 
also in agreement (Fig.|4]), though we do note that there are discrep- 
ancies in the axis ratios and external shear parameters, especially 
between the NICMOS and H-band models (see Section [57X1) . How- 
ever, the precision with which the lens parameters can be recovered 
is significantly higher for the Keck AO data, with the highest pre- 
cision obtained for the K'-Band. This clearly shows that the higher 
resolution provided by the AO imaging allows for a more precise 
lens modelling, when compared to the HST imaging, despite the 
lower signal-to-noise ratio (SNR). 

Finally, it is interesting to note that all of the models using 
the adaptive and grid-based Bayesian technique lead to a total mass 
centroid position and flattening that are consistent with those of the 
parametric mass modelling presented in Section 14.11 Also, all of 
the models have a total density profile that is very nearly isothermal 
(i.e., 7 = 2) out to the Einstein radius of the system. 



5 DISCUSSION 

We now discuss our results. First, we give a comparison of mod- 
elling the gravitational lens system using high angular resolution 
data from space- and ground-based telescopes. We then discuss the 
properties of the lensing galaxy and the magnification of the back- 
ground source that we infer from our model. Finally, we investigate 
the possibility of there being a second gravitationally lensed source 
for this system. 



5.1 Comparing HST and AO data sets for gravitational lens 
modelling 

The precision and accuracy of mass model constraints from ex- 
tended sources on strong gravitational lens galaxies is set by three 
conditions: (i) the number of independent resolution elements 
across the lensed images (and their multiplicity), (ii) the average 
SNR of the lensed-image surface brightness distribution and (iii) 
the level of surface brightness structure in the lensed source. While 
the first two conditions set the level of information contained in 
the lensed images and how well, in principle, the mass model of 
the lens galaxy can be constrained, the third condition determines 
the level of covariance (i.e., degeneracy) in the mass model, where 
more structured sources in general lead to a lower level of covari- 
ance between model parameters. One of the goals of this paper has 
been to illustrate points (i) and (ii) by comparing the results of high 
SNR ffiTF160W-band data with both lower resolution and pixel- 
sampling, and lower SNR Keck AO imaging data with higher reso- 
lution and sampling. We implicitly assume that the intrinsic source 
structure in the F160W-, H-, and K'-bands is very much correlated. 

Despite the fact that the results of the modelling of the three 
very different data sets are promisingly similar and the errors are 
small, we do note that there are differences. In particular, there 
are discrepancies between the axis ratio and position angle (PA) 
of the lens mass distribution, and the external shear strength and 
PA (see Fig. HI). Some tension between these quantities can be ex- 
pected, as they are strongly covariant (i.e., the external shear mim- 
ics the flattening of the lens potential/mass-distribution) and lack 
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Table 2. Individual parameters of the grid-based reconstruction lens models. The first row of a given model (specified by the first column) represents the 
maximum-likelihood model solution, while the second row gives the mean values for each parameter's posterior probability distribution. The third and fourth 
rows show, respectively, the 68 and 95 per cent confidence intervals of the means. The model parameters themselves are described as follows, b is the model 
lens strength (not the Einstein radius as defined for a SIS mass model) in arcseconds. is the position angle of the mass distribution of the lensing galaxy (in 
degrees east of north) while q is the axis ratio. 7 is the power-law slope of the mass profile. F and Fq are, respectively, the magnitude and position angle (in 
degrees east of north) of an external shear source. Additionally, the global Bayesian evidence of the model is presented in the final column. 



Model 


b 


6 


q 


7 


r 




Evidence 


Mk 


0.452 


-22.4 


0.853 


2.05 


0.014 


-u.i 






0.413 


-24.3 


0.846 


2.12 


0.013 


-72.7 


48806.4 




[0.410,0.415] 


[-24.5,-24.1] 


[0.845, 0.848] 


[2.11,2.12] 


[0.012,0.013] 


[-73.2,-72.2] 






[0.408,0.418] 


[-24.6, -23.9] 


[0.843, 0.849] 


[2.10,2.13] 


[0.011,0.014] 


[-74.8, -71.6] 




M h 


0.447 


-22.3 


0.853 


2.05 


0.019 


-78.5 






0.410 


-28.5 


0.786 


2.12 


0.029 


-89.1 


46280.0 




[0.397,0.426] 


[-29.8,-27.3] 


[0.770,0.801] 


[2.09,2.15] 


[0.025,0.032] 


[-97.9,-84.1] 






[0.360,0.461] 


[-31.7,-25.9] 


[0.761,0.811] 


[2.03,2.22] 


[0.022,0.038] 


[-101.0,-73.4] 




M H st 


0.439 


-23.0 


0.856 


2.07 


0.016 


-74.1 






0.424 


-25.0 


0.918 


2.09 


0.031 


-42.4 


7027.16 




[0.372,0.454] 


[-38.7,-16.8] 


[0.890,0.957] 


[2.04,2.20] 


[0.009, 0.044] 


[-53.8,-32.8] 






[0.355,0.520] 


[-39.8,-8.0] 


[0.777,0.988] 


[1.93,2.23] 


[0.006, 0.046] 


[-59.0,-6.6] 





of information combined with errors on the data can lead to bi- 
ases in the maximum-likelihood solutions and the posterior prob- 
ability distributions of individual parameters. This is especially 
true for B 193 8+666, where a nearly circular mass distribution and 
largely featureless Einstein ring conspire to add uncertainty to the 
lens model. Specifically, the circularly symmetric mass distribution 
(coupled with nearly coaxial foreground and background galaxies) 
makes it difficult to constrain the lensing galaxy's mass slope - 
increasing the covariances between other parameters - while the 
smooth light distribution lacks the contrast needed to differentiate 
between small-scale variations on the model. Thus, even moderate 
variations in parameter space can leave the lensed light distribution 
relatively unchanged, allowing significantly different models to fit 
the data equally well. 

Looking at the marginalized probability distributions in Fig. 
|4] it is obvious that there is a strong covariance between the lens 
strength (b) and mass slope (7), especially in the H-band and NIC- 
MOS models. This is why model discrepancies are the strongest 
between these data sets. However, we note that these degeneracies, 
and the parameter uncertainties in general, are much smaller from 
the higher-resolution (but lower SNR) Keck AO data. This implies 
that at least for some systems, ground-based AO data can "out- 
perform" space-based data in precision - and most likely also in ac- 
curacy - due to the better sampled lens image structure. Of course, 
this does not suggest that higher SNR data is unnecessary: of the 
two AO-based models, the K'-band data (with its relatively higher 
SNR) provides much tighter constraints on individual parameters 
than the H-band model, and indeed, even breaks the mass-slope 
degeneracy found in the H-band and NICMOS models. Rather, the 
comparison simply shows that, between SNR and image resolution, 
resolution seems to dominate any limitations on model precision. 

One has to be careful, though, when attempting to general- 
ize this conclusion to other lens systems, because the result de- 
pends critically on the complex interplay between SNR, spatial 
resolution, and source structure. It could very well be that lens 
systems with highly- structured sources (e.g., space-based B- or U- 
band data, where star formation could be strong and the sources 
could be more structured) could out-perform the higher resolution 
K'-band AO data. This interplay is currently under study. In cases 
where there is similar data, and a bright star is available, however, 



our results show that ground-based AO data can perform signifi- 
cantly better than their space-based HST data sets in constraining 
lens models. As we show in the companion paper, the higher reso- 
lution provided by AO imaging was also very power ful in the dis- 
cover y of a low-mass substructure in this lens system dVegetti et al.l 
120121) . 



5.2 The foreground lensing galaxy 

The B 1938+666 lensing galaxy is clearly seen in the optical and 
NIR data sets (Fig. Q]), and its surface brightness distribution is 
well-modelled by a single elliptical Sersic profile (see Section 
14.11 ). Excluding the HST V-band photometry, we find the galaxy 
to have I— H and H— K' colours of 1.9 and 0.4 magnitudes, re- 
spectively, after Galactic reddening corre ctions. We model these 
colours with the lBruzual & Charlotl (12003b stellar population code 
and, assuming a Solar metallicity and no dust, find consistency 
with an old (> 4 Gyr) stellar population. The foreground galaxy's 
isophotal regularity, colours, and a bsorption line spectrum mea- 
sured by iTonrv & Kochanekl d2000h suggest that the lens has an 
early- type morphology; this is consistent with the previous analy- 
sis of iKing et alJ ll998h based only on the NICMOS F160W-band 
imaging. 

We use the lBruzual & C harlot (20031) 4 Gyr stellar population 
model to compute the V-band absolute magnitude and to estimate 
the stellar mass within the Einstein radius. The luminosity is found 
to be Lv = 3 x 10 10 L© and the stellar template has a V-band 
stellar mass-to-light ratio of 1.54, implying a stellar mass within 
the Einstein radius of M* = 4.7 x 10 10 M© if a Chabrier ini- 
tial mass function (IMF) is assumed. We can compare this to the 
total mass within the Einstein radius from lensing, and we find a 
dark-matter mass fraction of Mdark/Mi en s = 0.55 within the Ein- 
stein radius. If we instead use a Salpeter IMF to describe the stel- 
lar population, the stellar mass within the Einstein radius becomes 
M* = 8.2 x 10 10 M©, leading to a dark-matter mass fraction of 

M dark /Miens = 0.2. 

Ideally, we would like to compute the stellar- to- total mass 
fraction within an aperture physically associated with the lensing 
galaxy. Therefore, we use two methods to estimate the effective 
radius from the more robustly measured lensing data. First, we em- 
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Figure 3. The gravitational lens mass model of the B 1938+666 Einstein ring, using the adaptive grid-based method of Vegetti & Koopmans < 2009a|). The 
results are for the three independent data sets: (a) the AO K'-band, (b) the AO H-band and (c) the HST NICMOS H-band. For each plot, we show the measured 
surface brightness distribution of the Einstein ring (top left panel), the best-fit smooth lens model reconstruction (top right panel), the residual image (bottom 
left panel) and the reconstructed unlensed image of the background source galaxy (bottom right panel). We note that the source reconstruction grids are not 
registered. Therefore, taking into account registration offsets, as well as small offsets due to PSF convolution, we find that the source positions are consistent 
with being coincident in all data sets. 



ploy the relationship be tween the power-law density slope and ef- 
fective radius found bv lAugeretaHd2010h and co rrect for early- 
type galaxy growth rates (e.g., New man et al There are sev- 
eral significant caveats, including evoluti on of the power- law slope 
fe.g jRuff et alJ201ll : lBolton et al.l2012h and large uncertainties on 
these relationships; nevertheless, we find that the effective radius is 
consistent with being the same size as the Einstein radius. Moti- 



vated by this, we check if the assumption that the effective radius 
of 0.45 arcsec (i.e., equivalent to the Einstein radius) is consistent 
with measurements of the fundamental plane. We use the veloc- 
ity dispersion inferred from the lensing model as a proxy for the 
stellar velocity dispersion and apply a p assive evolution correction 
to the luminosity as determined by the iBruzual & Char iot (2003) 
4 Gyr template. We find that, again correcting for size evolution, 
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Figure 4. The marginalized posterior probability distributions of the lens model parameters, as measured by a nested sampling analysis. Each grid shows a 
different marginalization: the 2-dimensional contours represent the distribution between two model parameters (specified by the row and column), while the 
1 -dimensional histograms at the top of each column represent the distribution of a single parameter. The AO K ; -band, AO H-band, and NICMOS data are 
represented by red, green, and blue contours, respectively. The maximum-likelihood model values for each parameter, and their uncertainties, are presented in 
Table[2] 



our assumption that the effective radius and Einstein radius are co- 
in cident yields consis tency with the fundamental plane presented 
bv lAuser et aDfeOlOh . 



5.3 The magnification of the NIR light from the background 
galaxy 

The gravitational lens models also provide an estimate of the source 
surface brightness distribution, after correcting for the lensing ef- 
fect. From the pixelated reconstruction, we find that the source is 
composed of two components at NIR wavelengths (in agreement 
with the parametrized source distribution discussed in Section flTTV 
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One component is a high surface brightness region with a projected 
size of ~ 0.8 kpc, the second component is more extended, up to 
~1.6 kpc in projected size, and has a lower surface brightness. We 
defer any detailed discussion about the nature of the NIR source 
reconstruction - especially with respect to the AGN and molecular 
gas components traced by the emission at radio wavelengths - to a 
follow-up paper. 

To estimate the total magnification of the NIR source from the 
gravitational lens models, we compare the model emission found 
in the lensing plane to that found in the source plane. We use the 
HST NICMOS data set for this calculation because these data have 
the most robust measurement of the extended light distribution of 
the source galaxy; the AO data sets will have some light artificially 
pushed into an extended envelope (thus leading to a greater chance 
of confusion between source and lens galaxy light) due to the fact 
that AO Strehl ratios are less than unity. Additionally, the brightness 
of the sky background at NIR wavelengths can hide faint, extended 
emission at the outer edge of the source galaxy, leading to a biased 
magnification estimate. 

Overall, we find that the total magnification factor of the NIR 
emission is ^ 13. We note that a magnification of 176 has been pre- 
viou sly reported for B 193 8+66 6 from a simple point- source model 
(see iBarvainis & IvT son 2002 for some details). This larger mag- 
nification was used by iRiechersI d201 lh for their analysis of the 
molecular gas properties of B 1938+666, and they found that the CO 
(3-2) line intensity for B 1938+666, and hence the molecular gas 
mass, was about an order of magnitude lower than for other quasars 
at a similar epoch. This large discrepancy is almost certainly due to 
the magnification of the gas emission being over-estimated. How- 
ever, although we find a smaller magnification for the NIR emission 
region for this source, some caution should be taken when applying 
this value to other wavebands, particularly if the emission has a dif- 
ferent size, or position with respect to the lensing caustics produced 
by the lens; only spatially resolved imaging and a good lens model 
can give a robust estimate of the magnification for any particular 
emission region of a gravitationally lensed source. 



5.4 Luminous Substructure 



In lVegettietaD <2012l) we presented the detection of a low-mass 
substructure in the B 1938+666 system from the distorting effect 
the substructure has on the lensed arc. We do not find any ev- 
idence of luminous, and probably more massive, substructure in 
the residual i mages for the system, which is consistent with the 
IVegetti etaD j2012h analysis. To determine an upper limit on the 
brightness of any luminous substructure, we add a series of sim- 
ulated point sources between 24.5 < m K t < 28 to an image of 
the lens system from which both the lens galaxy and the Einstein 
ring emission have been subtracted. For our point source model, 
we choose the empirical PSF star used in the lens modelling, as 
this would accurately represent how a point source object would 
appear in the image. One example simulation is presented in Fig. [5] 
We choose a point source model over other models - such as Gaus- 
sian or Sersic profiles - because we have no inherent knowledge 
about the morphologies of these objects, and because the typical 
small size of satellites relative to a parent galaxy, coupled with the 
high redshift of the B 1938+666 lensing galaxy, suggests that these 
objects would be likely unresolved in our data. Additionally, since a 
point source object has a higher surface brightness than an extended 
object of the same magnitude, a point- source limiting magnitude 
represents a robust limit to the lowest luminosity that can possibly 







Figure 5. Residual image of the B 1938+666 system with simulated point 
source objects added, which we use to estimate a detection limit for lumi- 
nous substructure. The point sources detected by SExtractor are represented 
by green, solid-line circles, while those that are missed are represented by 
red, dashed-line circles. There is a distinct difference in noise level between 
the regions inside and outside of the Einstein ring, suggesting that the sub- 
structure limiting magnitude should be brighter closer to the lens. 



be observed, placing a hard limit on the luminous properties of any 
substructure that is independent of morphology. 

We generate a total of 1000 simulated images, randomizing 
both the magni tudes and positions of th e point sources, and then 
use SExtractor dBertin & ArnoutslH996l) to detect the added satel- 
lites. We consider a point source to be detected if SExtractor lo- 
cates an object within 5 pixels and 1.5 magnitudes of the actual 
position and magnitude of the simulated object. Compiling the re- 
sults into magnitude bins, we find that we are able to recover the 
majority (i.e. > 50 per cent) of the simulated objects down to an 
apparent magnitude of K' = 26.1, which we take to be our de- 
tection limit (see Fig. [6]). Converting this value to the rest-frame 
V-band corresponds to an object that has an absolute magnitude 
of My = —16.2, which is fainter than the Magellanic clouds 
(M^lmc = —18.5, Mv,smc = —17.1), but approximately three 
times brighter than t he Sagittarius dwarf satellite (Mv,s g r = — 15; 
iTollerudet ai1l2008l) . 



5.5 A second lensed source? 

There are two faint arclet structures that can be seen just outside 
the Einstein ring on the eastern side of the K' image (Fig. Q] bot- 
tom middle panel). To enhance the signal of these arclets, we apply 
a Gaussian filter (a =6.5 pixels) to smooth the data, followed by 
a Laplacian filter to increase the contrast around bright flux peaks. 
Fig. [7] shows a high-contrast image of the original data and an im- 
age showing the results of applying Gaussian and Laplacian filters. 
While the two arcs are barely detectable above the noise in the un- 
processed data, the arcs are clearly visible in the filtered image. 
As the arclets are so faint in the original K'-band data, we do not 
formally consider them in our lens model. However, taking into ac- 
count the enhanced image, it is possible - given their positions and 
orientations relative to the lensing galaxy - that these arclets are 
actually multiple images of a second background source. 
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Figure 7. Left: K'-band image of B 1938+666 with the two faint gravitational arc-like structures highlighted (A and B). Right: same image after applying a 
Gaussian filter to smooth the image, and then a Laplacian filter to highlight the faint features. Both potential gravitational arcs can clearly be seen. 




24.5 25.0 25.5 26.0 26.5 27.0 27.5 28.0 
Magnitude 



Figure 6. Fraction of simulated point sources recovered by SExtractor, 
binned as a function of magnitude. We are able to recover over 50 per cent 
of sources in each magnitude bin, up to and including the m K i = 26.1 
bin (represented by the dashed line). We therefore treat this as our limiting 
magnitude. 



There are three possible scenarios: the arclets are not being 
lensed; the arclets are images of a strongly-lensed source that is as- 
sociated with the primary ring (e.g., an additional component of the 
galaxy being lensed into the ring, a nearby galaxy, or more images 
of the ring - although this would lead to an extremely peculiar im- 
age configuration that we do not expect to arise from well-behaved 
lensing mass distributions); or the arclets are images of a second 
source being strongly lensed but at a completely different redshift 
than the Einstein ring galaxy. We can investigate these scenarios by 
using the lens models that have been fitted to the Einstein ring. We 
first use our best lens model to cast the emission from the Einstein 
ring and the two arclets back to the source plane. As expected, the 
arclets do not fall on the Einstein ring source. However, the two 
arclets also are not mapped to a coincident location on the source 
plane, indicating that if they are at the same redshift as the Einstein 
ring galaxy they are not multiple images of the same source. 

The possibility that the arclets correspond to images of a 



strongly-lensed source at a different redshift is intriguing; such 
double-plane lens systems provide very tight constraints to the 
mass-density slope of the len s and can give usefu l constraints on 
cosmological parameters (see iGavazzi et al.l 12008 for details). To 
investigate this possibility we must rescale the lens strength of our 
mass model, where the scaling is directly proportional to the ratio 
of angular diameter distances between the lens and the two sources. 
We find that the two arclets can be mapped to approximately the 
same location in the source plane if the lens strength is increased 
by a factor of ~ 1.85. Such a large scaling factor requires an unreal- 
istically high redshift for the source, although two caveats exist that 
could mitigate this concern: the galaxy being lensed into the Ein- 
stein ring can contribute to the lensing, and the mass profile could 
deviate from the central power law that was found from the Einstein 
ring fit. 

In any case, some estimate of the redshifts of the two poten- 
tial gravitational arcs will need to be obtained to determine their 
nature. This will likely require extensive followup with deep, high- 
resolution spectroscopy. While this is difficult with current tele- 
scopes (especially because arclet B is so faint), future instruments 
such as the James Webb Space Telescope (JWST), the Thirty Me- 
ter Telescope (TMT), or the European-Extremely Large Telescope 
(E-ELT) should provide the resolution and light-gathering power 
necessary to achieve this goal. 



6 CONCLUSIONS 

We have presented a new mass model for the gravitational lens sys- 
tem B 1938+666, using a grid-based Bayesian reconstruction tech- 
nique on high resolution ground-based AO and space-based HST 
data, as part of the Strong-lensing at High Angular Resolution Pro- 
gram (SHARP). We find that the smooth component of the lens- 
ing galaxy's mass profile is well-fitted by a (nearly isothermal) 
power-law distribution, while its light profile (and that of the source 
galaxy) can be described by Sersic components. A more in-depth 
analysis of the mass model, characterizing the amount of substruc- 
ture prese nt in the system, is presented in a companion paper to 
this work (Vegetti et al. 2Ql3). The model is consistent over three 
independent NIR data sets - due largely to the constraints provided 
by a bright Einstein ring - and agrees well with previously reported 
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results. When compared to models derived from traditional ground- 
based imaging, though, the high-resolution models are significantly 
more precise. 

The relative improvement in precision varies from data set to 
data set. Thus, by generating both AO- and ffiT-based models, 
we have provided a quantitative comparison between instruments. 
Overall, we find that the uncertainties on the model parameters de- 
rived from the AO data sets are smaller than those measured from 
the HST, by as much as an order of magnitude. This suggests that 
AO data are better at constraining lens models than equivalent HST 
data (at least in cases where the lens and source galaxies are red and 
have reasonably smooth light profiles), and therefore, that AO ob- 
servations of lens systems can lead to a better description of those 
systems' mass distributions. However, confirmation of this result 
will require a much larger data set and will be explored in a future 
paper (SHARP II; Fassnacht et al., in prep). 
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